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*  Approximate  analytical  expressions  for  electron  distribution  functions  in  III-V  semi¬ 
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two  analytic  methods  were  used  and  compared  with  an  earlier  iteration  method  and  a  new 
perturbation  method.  All  gave  similar  distribution  functions  and  drift  mobility.  For 
high  electric  fields  the  electron  distribution  function  was  found  within  the  diffusion 
approximetion  by  an  expansion  Inversely  in  the  electric  field.  For  A1q  25In0  75  the 
effects  of  p-wave  scattering  and  nonparabolic  energy  bands  were  not  fotma'to  oe  the  cause 
of  the  negative  differential  resistance.  Alloy  scattering  in  ternary  semiconductors  was 
treated  in  the  coherent  potential  approximation.  The  validity  of  the  diffusion  and  maxi¬ 
mum  anisotropy  approximation  has  been  asseesed  by  an  expansion  inversely  in  the  electric 
field  of  the  polar  optical  collision  term  in  high  electric  fields .a 
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I.  OBJECTIVES  OF  RESEARCH 


The  objective  of  this  research  was  to  derive  approximate  analytic  expressions 
for  electron  distribution  functions  in  III-V  semiconductors  in  order  to  calculate 
transport  properties*  The  approach  was  viewed  as  an  alternative  to  Monte  Carlo  and 
drifted  Maxwellian  calcuatlons.  The  use  of  analytic  expressions  for  the  electron 
distribution  function  was  intended  to  reduce  the  cost  of  numerical  calcuatlons  and 
to  facilitate  parameter  studies  of  transport  coefficients 9  in  particular  the  temper¬ 
ature  dependence  of  the  drift  mobility*  The  objective  during  the  two  years  of  this 
contract  was  to  focus  on  the  high  electric  field  regime  for  GaAs,  InP  and  ternary 
semiconductors  such  as  AIq  25  Iq  75  As.  For  ternary  III-V  semiconductors  nonparabolic 
energy  bands  were  to  be  included  in  the  central  valley  for  polar  optical  scattering 
in  order  to  understand  the  negative  differential  resistance  predicted  for  electrons 
in  the  central  valley  alone.  Alloy  scattering  was  also  to  be  considered*  An  at¬ 
tempt  to  Include  Intervalley  scattering  for  GaAs  was  to  be  made*  All  new  results 
were  to  be  compared  with  previous  numerical  and  experimental  results  where  possible. 
II.  SUMMARY  OF  RESEARCH  RESULTS 

Several  related  research  topics  have  been  explored.  The  main  thrust  of  the 

research  has  been  the  development  of  theoretical  methods  for  calculating  electron 

distribution  functions  and  the  drift  mobility  in  III-V  semiconductors.  The  contract 

has  supported  in  addition  to  the  principal  investigator  one  graduate  student 9  Sik 

Han  Soh,  who  completed  his  master 9 s  thesis  in  August 9  1984.  He  is  now  employed  as 

a  project  engineer  working  with  III-V  semiconductors  at  Hewlett-Packard  in  Singapore. 

His  thesis1  was  entitled  "Electron  Distribution  Functions  and  Mobility  in  III-V  Seml- 

2 

conductor  Compounds  for  High  and  Low  Electric  Fields".  A  paper  on  part  of  this  work 
appears  in  The  Physics  of  Subnicron  Structures  which  has  been  published  by  Plenum. 

A  second  paper  is  planned  on  the  effect  of  p-wave  scattering  and  nonparabolic  energy 
bands  on  the  negative  differential  resistance  arising  in  the  central  valley  of  the 
conduction  bend  for  certain  ternary  III-V  semiconductors.  Also  a  short  collaboration 


with  Professor  Peter  Fedders  of  the  Physics  Department  at  Washington  University 

and  Professor  Charles  Myles  of  the  Physics  Department  at  Texas  Tech  University  on 

the  coherent  potential  approximation  applied  to  alloy  scattering  In  ternary  III-V 

3 

semiconductors  resulted  In  a  talk  .  For  comparable  percentages  of  disorder  the 
scattering  due  to  off-diagonal  disorder  was  found  to  dominate  In  alloy  scattering 
in  a  two-band  model  for  III-V  semiconductors.  This  project  was  not  pursued  further 

by  the  principal  Investigator  because  a  subsequent  paper  by  Professors  Dan  Rode 

4  * 

and  Peter  Fedders  showed  that  for  most  ternary  III-V  semiconductors  of  Interest 
in  device  applications  the  Born  approximation  could  be  justified  such  that  the 
coherent  potential  approximation  Is  not  really  needed.  A  resolution  of  some 
of  the  differences  in  the  diffusion  approximation^  and  the  maximum  anisotropy 
approximation^ f 7  used  to  simplify  the  electron  distribution  function  at  high  d.c. 
electric  fields  has  been  made  by  means  of  an  asymptotic  expansion  inversely  In 
the  electric  field  of  the  steady-state  Boltzmann  equation  for  polar  optical  scatter¬ 
ing.  A  better  understanding  of  the  applicability  of  these  approximations  Is  now 

g 

feasible  .  Preliminary  calculations  on  the  extension  of  the  calculation  of  the 
electron  distribution  function  by  approximate  analytic  expressions  for  intervalley 
scattering  was  begun.  That  effort  was  not  continued  because  of  an  uncertainty  in 
the  basis  of  the  validity  of  the  diffusion  and  maximum  anisotropy  approximations 
and  because  the  principal  Investigator  ran  out  of  time  in  one  year.  The  drift 
mobility  for  InP  was  calculated  but  not  investigated  further  as  the  time-of-fllght 
data  was  not  available. 

The  master's  thesis  of  Mr.  Soh  included  several  methods  for  determining 
electron  distribution  functions  in  low  electric  fields  and  the  calculation  of  the 

Q 

electron  drift  mobility  for  polar  optical  scattering.  The  methods  were:  iteration, 
perturbation,  and  two  analytic  methods.  The  iteration  method  was  calculated  nu¬ 
merically  and  compared  with  previous  numerical  results  that  had  been  compared  with 
experimental  data.  The  perturbation  method  was  numerical  but  the  two  analytic 


methods  produced  approximate  analytic  expressions  for  the  electron  distribution 

function  from  which  the  drift  mobility  was  calculated.  The  results  from  all  methods 

were  similar.  For  high  electric  fields  the  expansion  of  the  steady-state  Boltzmann 

equation  inversely  in  the  electric  field  within  the  diffusion  approximation  for  the 

2 

electron  distribution  function  was  reviewed  .  The  electron  distribution  function 
was  compared  with  that  derived  by  Stratton^  for  parabolic  energy  boards  and  agreement 
found  for  an  expansion  in  l/£  of  his  result.  Neither  distribution  function  is 
normalizable  and  intervalley  scattering  needs  to  be  included  for  GaAs;  for  certain 
ternary  III-V  semiconductors  intervalley  scattering  may  be  omitted  in  a  certain  range. 
For  high  electric  fields  p-wave  scattering  and  nonparabolic  energy  bands  were  in¬ 
cluded  in  the  polar  optical  scattering  term.  For  GaAs  and  even  A1q  ^5  75  A® 

the  effects  of  p-wave  scattering  and  nonparabolic  energy  bands  were  not  found  to 
alter  the  electron  distribution  significantly.  Consequently  the  negative  differential 
resistance  for  ternary  III-V  semiconductors  does  not  appear  to  be  caused  by  the  non- 
parabolic  character  of  the  central  valley  of  the  conduction  band.  Another  supporting 
argument  for  our  conclusion  can  be  found  in  the  Monte  Carlo  calculations^  of  the 

electron  drift  velocity-electric  field  curve  where  p-wave  scattering  is  negligible. 

12 

As  Rode  has  pointed  out  nonparabolic  energy  bands  exist  only  if  there  are  p-wave 
matrix  elements  as  well  as  s-wave  matrix  elements.  Hence t  one  cannot  Invoke  non- 
parabolic  energy  bonds  at  the  same  time  p-wave  scattering  is  dropped. 

To  compare  the  diffusion  and  maximum  anisotropy  approximations  an  expansion 
inversely  in  the  electric  field  has  been  made  of  the  steady-state  Boltzmann  equation 
for  a  d.c.  electric  field  and  polar  optical  scattering.  This  model  is  expected  to 
be  reasonable  for  ternary  semiconductors  if  alloy  scattering  is  added;  for  GaAs 
intervalley  scattering  should  be  added.  For  E0/R  small  where  Is  the  critical 
electric  field  characterised  by  the  strength  of  the  coupling  to  the  polar  modes, 
the  symmetric  part  of  the  distribution  function  for  electrons  in  the  central  valley 
is  a  constant.  The  asymmetric  part  of  the  electron  distribution  function  can  be 


represented  as  an  expansion  in  odd  order  Legendre  polynomials  whose  argument  is 
the  cosine  of  the  angle  between  the  electric  field  and  wave  vector.  The  coefficients 
of  the  odd  order  Legendre  polynomials  for  high  electric  fields  and  electron  energies 
high  compared  to  an  optical  phonon  energy  obey  a  modified  maximum  anisotropy  relation 
where  no  even  order  coefficients  for  n  >  0  exist.  The  diffusion  approximation 
does  not  give  the  entire  electron  distribution  function  correctly  but  for  high 
electric  fields  it  produces  a  functional  from  which  the  electron  density  and  drift 
velocity  for  electrons  in  the  central  valley  can  be  correctly  calculated  for  polar 
optical  scattering. 
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